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. . The framework for considering the astronomical and cosmological observations in the context 

04 , of scalar-tensor quintessence in which the quintessence field also accounts for a time dependence of 

(«_^ ■ the gravitational constant is developed. The constraints arising from nucleosynthesis, the variation 

^-^ ' of the constant, and the post-Newtonian measurements are taken into account. A simple model of 
supernovae is presented in order to extract the dependence of their light curves with the gravita- 

51^ ' tional constant; this implies a correction when fitting the luminosity distance. The properties of 

^ ' perturbations as well as CMB anisotropies are also investigated. 
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(^ , I. INTRODUCTION 

^-^ ' The combination of recent astrophysical and cosmological observations (among which are the luminosity distance 

versus redshift relation up to z ~ 1 from type la supernovae Q], the cosmic microwave background temperature 
anisotropies and gravitational lensing 0]) seems to indicate that the universe is accelerating and thus that about 
i-^ ■ 70% of the energy density of the universe is made of a matter with a negative pressure (i.e., having an equation 
hr^i of state uj = P/p < 0). This raises the natural question of the physical nature of this matter component. Indeed, 
Q ■ a solution would be to have a cosmological constant (for which uj = —1) but one will then have to face the well 
H ' known cosmological constant problem El, i.e., the fact that the value of this cosmological constant inferred from the 
t^ ' cosmological observation is extremely small — about 120 order of magnitude — compared with the energy scales 



^ 



of high energy physics (Planck, grand unified theory, strong and even electroweak scales). Another solution is to 
argue that there exists a (yet unknown) mechanism which makes the cosmological constant strictly vanish and to 
find another matter candidate able to explain the cosmological observations. Indeed it assumes that the cosmological 
constant problem is somehow solved and replaces it by a dark energy problem, 
d ' This latter route has gained a lot of enthusiasm in the past years and many candidates have been proposed (for 

recent reviews on see, e.g., Refs. [HH]). Among all these proposals, quintessence |^-|^ seems to be the most promising 
mechanism. In these models, a scalar field is rolling down a potential decreasing to zero at infinity (often referred to 
as a runaway potential) hence acting as a fluid with an equation of state varying in the range — 1 < w < 1 if the field 
is minimally coupled. Runaway potentials such as the exponential potential and inverse power law potentials 

n^) - -^^^ (1) 

with a > and M a mass scale, can be found in models where supersymmetry is dynamically broken |1G|] and in 
which flat directions are lifted by nonperturbativc effects. 

As clearly explained in Ref. [|ll[, all the models for this dark energy have to (i) show that they do not contain in 
a disguise way a cosmological constantlike fine-tuning (the fine-tuning problem), (ii) explain why this kind of matter 
starts to dominate today (the coincidence problem), (iii) give an equation of state compatible with the observational 
data (the equation of state problem) and (iv) arise from some high energy physics mechanisms (the model buil ding 
problem) . Quintessence models mainly solved the fine-tuning problem because of the existence of tracking solutions 112] 
(first studied in Ref. [Rjq] ) which are scaling attractor solutions of the field equations and allows the initial conditions 
for the scalar field to vary by about 150 orders of magnitude. The second tuning (related to the coincidence) concerns 



the mass scale M that has to be determined by the requirement that about 70% of the energy density of the universe 
is in the quintessence field. As shown in Ref. [|l3| for the case of the inverse power law potential, this mass scale is 
comparable to other scales from high energy physics and the tuning on this mass scale is mild provided the exponent 
a is not too small (a must be bigger than 4 so that M > 1 TeV). The equation of state depends on the shape of the 
potential and it is hoped that it will be soon determined by, e.g., a weak lensing experiment p4| . For instance, it has 
been shown that an exponential potential cannot lead to an accelerating universe ||7|,g| and that the equation of state 
for an inverse power law potential mainly depends on the slope a. As explained above, the model building is also 
addressed in the framework of supersymmetry 

Quintessence scenarios, however, have some important problems. The requirement of slow roll (mandatory to have 
a negative pressure) and the fact that the quintessence field dominates today imply that (i) it is very light |l^] 
(roughly of the order ~ 10"'^^ eV and it should induce violation of the equivalence principle and time variation of the 
gravitational constant) and that (ii) the vacuum expectation value of the quintessence field today is of order of the 
Planck mass. This latter problem led Brax and Martin [|ll|,f6J to propose that a supergravity correction had to be 
taken into account leading to the so-called supergravity (SUGRA) quintessence potential 
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which shares the same properties as the inverse power law potential at early time but which stabilizes the quintessence 
field accounting for a better agreement of the equation of state. 

Note that the two main features of a quintessence potential is that it must be steep enough for the field to be in 
a kinetic regime for a large set of initial conditions, hence redshifting faster than radiation and being subdominant 
at nucleosynthesis, and then to reach a slow roll regime to mimic a cosmological constant. This latter regime will 
always ultimately take place and the parameters of the potential have to be tuned so that it happens around today. 
The simplest waya to implement this idea are to use inverse power law potentials and exponential potentials which 
are one parameter potentials. Another solution (but involving more parameters) is to consider potentials with a local 
minimum as was first proposed by Wetterich B . 

An underlying motivation to replace the cosmological constant by a time dependent scalar field probably lies 
in string models in which any dimensionful parameter is expressed in terms of the fundamental string mass scale 
and the vacuum expectation value of a scalar field. For instance, string theory has revived the consideration of 
gravitational-strength scalar fields [|l7| such as Kaluza-Klein moduli or the dilaton appearing in the low energy limit 
of the gravitational sector leading to scalar-tensor theories of gravity. As explained above, the quintessence field is 
expected to be very light and this points toward scalar-tensor theories of gravity in which a light (or massless) scalar 
field can be present in the gravitational sector without being phenomenologically disastrous. These arguments lead 
to consider quintessence models in the framework of scalar-tensor gravity. Indeed, the dilaton and the quintessence 
field can be two different scalar fields (as considered, e.g., in Ref. |l^ in the particular case of Brans-Dicke theory) 
or the same scalar field. The latter subclass involving a single scalar field (the quintessence field is also the dilaton) 
dictating the time variation of both the gravitational constant and the cosmological constant is attractive; it involves 
less free functions and has been focused on much in recent years. The study of these quintessence models, referred to 
as nonminimal quintessence |I^ , coupled quintessence [^ , extended quintessence |^ or generalized quintessence p^ ] 
was mainly motivated by the fact that tracking solutions have been shown to exist for nonminimally coupled scalar 
fields |l|,|2|l. 

Scalar-tensor theories are the most natural extensions of general relativity, in particular they contain local Lorentz 
invariance, constancy of nongravitational constants and respect the weak equivalence principle. The most general 
action for these theories [p3[ for the matter and gravity is given, in the Jordan frame, by 
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where £,-„at is the Lagrangian of ordinary matter (such as radiation and pressureless matter), o is the determinant of 
the metric g^j^, R the Ricci scalar, and V^ is a potential to be discussed below. The action (|3|) depends a priori on 
three arbitrary functions F, uj, and V. If F{R, cj)) is not a trivial function of R then one has an additional (massive) 
scalar degree of freedom |2J] so that the requirement that there is only one scalar partner to the graviton implies 



^(i?,0) = — F(0)i?, 



(4) 



where k = 8ttG = A/pj, G being the bare Newton constant. F{(p) is a dimensionless function which needs to be 
positive to ensure that the graviton carries positive energy p5[. Now, by a redefinition of (f> we can always choose 



either F{(j)) or u!{(j))/(j) to be unity so that we are left with only two independent functions of the scalar field one being 
the potential V . Choosing lu{(J>)/4> = 1, the effective (time dependent) constant is related to G by 

GM - ^. (5) 



As was explained above, authors first considered models with a nonminimally coupled scalar field [|19| , p6H 2q|, i.e., 
in which the function F is decomposed as 

Fi,p) = l + 2^Kf{(t>). (6) 

Indeed, as long as we have not fixed the choice of the function /, such a decomposition is completely general and 
always possible. If the normalization of / is chosen in a way that /(Mpi) ~ ©(MpJ today then ^ gives an order of 
magnitude of the deviation with respect to general relativity today. 

Due to the time variation of the gravitational constant, the strength of the coupling ^ can be constrained once 
the function /(</)) has be chosen. Chiba |E9[ gave the constraints arising from the post-Newtonian (PN) parameters 
and the time variation of G for / = (jr/2. In that case the deviation of the scalar-tensor theory from general 
relativity was fixed but, as pointed out by Bartolo and Pietroni |3C[| in scalar-tensor quintessence models a "double 
attractor mechanism" can happen, namely, of the scalar-tensor theory towards general relativity (through the Damour- 
Nordtvedt mechanism [|l|) and of the quintessence field toward its tracking solution hence allowing for large deviation 
from general relativity in the early universe. Note that in these models the dilaton is not completely stabilized and 
is slow rolling in its runaway potential. 

Many works have then studied the cosmological implications of these models, starting from the study of the 
perturbations in Brans-Dicke theory [g2[ and for a nonminimally coupled scalar field |^,Q the computation of 
cosmic microwave background (CMB) anisotropics [gl],^, the properties of the nucleosynthesis |p0| , p6| . But among 
all, one of the very interesting results concerns the possibility to rule out some of these models p3] : it was shown that 
because of the positivity of the energy of the graviton, one can mimic a model in which the gravity is described by 
general relativity with a cosmological constant by a scalar-tensor theory with V = (or when certain relation between 
V and F are set) only in a certain range of redshift, hence offering a powerful test on a large class of scalar-tensor 
models. 

In this article, we try to extract the observational constraints on these quintessence models in scalar-tensor gravity. 
We first recall in Sec. the background properties: we study the constraints arising from the bounds on the post- 
Newtonian parameters, the time variation of the gravitational constant, of nucleosynthesis and of the positivity of the 
energy of the graviton. We then turn to the information that can be extracted for the type la supernovae (Sec. [II B| ). 
This lead us to describe a simple model of SN la in order to extract the dependence of the light curve with the 
gravitational constant and give us a modified magnitude vs redshift relation that has to be used when extracting 
the luminosit y dis tance vs redshift relation. We finish by a computation of the angular diameter distance vs redshift 
relation (Sec. II C ). We then turn to the property of the perturbations and generalize the attractor property (Sec. |l|) 



found in Ref. [|37| to the scalar-tensor quintessence and then discuss the CMB angular spectrum (Sec. IV). All the 
technical details are gathered in the appendices and we apply these results all along the article to two examples, one 
being a nonminimally coupled scalar field and the second a scalar-tensor theory with conformal coupling. 

II. BACKGROUND PROPERTIES 

We consider a Friedmann-Lemaitre universe with the line element 

ds2 = dt^ - a^{t)jijdx'dx^ = g^^dxf'dx'', (7) 

where a is the scale factor, t the cosmic time and jij the metric on constant time hypersurfaces. Greek indices run 
from to 3 and Latin indices from 1 to 3. As detailed in Appendix ^ the Friedmann equation in presence of a 
nonminimally coupled scalar field takes the form 

H'' + K= ^a2 (p„„, + p,,e + P,) , (8) 

where we have introduced the conformal Hubble parameter Ti. = a/a, and where an overdot denotes a derivative with 
respect to the conformal time {dt — adry), and K is the curvature index. On the right-hand side, p„^^t designs the 
matter energy density, and Pmc and p^ , the energy densities of the scalar field and respectively defined respectively in 



Eqs. ( [B^ ) and ( B10| ). The quantity k^^ = k/F{(J)) acts as an effective Newton constant and depends also on (j) (and 



hence, varies with time). We define the density parameter of any component / by 

where one has to use k^{{ and not k, since we have only access to a measure of the effective Newton constant and not 
of its bare value. 

A. How is F{(j>) constrained? 

In this section, we review the different constraints on F{(j)) arising from the requirement of the positivity of the 
energy of the graviton, the bounds on the post-Newtonian parameters and of the time variation of the gravitational 
constant and from nucleosynthesis. 



1. Positivity of the energy of the graviton: Diagonalizing the action (M) by the conformal transformation |25 



2 - 2 



one can show that the graviton is the perturbation of g^i/ and its scalar partner evolving in the potential Ujcj))- 
It can be shown that a scalar-tensor theory in fact well defined only if the transformation (nfl) is possible p3,p5[ . 
This is the case if F(0) is positive. When integrating our equations this constraint has to be checked. In 
practice, since one has k^u cx F~^ cx Ti^ -f K, this constraint is always satisfied when the universe is spatially 
flat. 

2. Post-Newtonian constraints : With the form (||), the standard post-Newtonian parameters are given by ||2^ , P5|| 
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where a dash denotes a derivative with respect to the field (j) and where a subscript means that the function 
is evaluated today. The function a is defined by 

,~, dln^ , , 

a((/.) = — ^. (13) 



Current constraints (see, e.g., Ref. |38| for a recent review of the measurements) give 

|7pn - 1| < 2 X 10~^ |/?PN - 1| < 6 X lO""*. (14) 

This implies that 
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o-mKFJ^o)<10-\ (15) 

and, as explained in Rcf. |22| , the second bound cannot be used to constraint dao/di/). 

3. Time variation of G : The effective constant Goff deduced from Eq. (0) is not the gravitational constant that 
would be measured in a Cavendish-Michel-type experiment, i.e., it is not the effective Newton constant. This 
constant GN^ff, entering in the force between two masses, is given by 

G,^^^ = ^'J&LtA(^^GAUl^al), (16) 

Fo2M^^Fo + 3{fX 



in which one has two contributions, namely the exchange of a graviton and of a scalar |25 
straints l39| on the variation of the Newton constant imply 
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4. Nucleosynthesis: Nucleosynthesis bounds have two origins. Roughly speaking, we have to require that (i) the 
matter contents dominating the Friedmann equation behaves as radiation at nucleosynthesis and that (ii) the 
effective number of degrees of freedom of the relativistic particles, g* say, does not vary from more than 20% 
than its expected value g^, — 10.75 at this epoch. 

In the case of quintessence with an exponential potential, these bounds were used to show that the quintessence 
field can not close the universe Gffl. For general inverse power law and SUGRA potentials, the first constraint 
was shown |13| to imply bounds on the initial condition of the scalar field at the end of inflation. Here, we want 
to estimate the bounds arising from the second requirement and we assume that the contribution of the scalar 
field is subdominant with respect to the radiation. The energy density of the radiation is given by 
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where T is the temperature. Assuming that pq ^ p^^^ [pq being defined by Eq. ( Bl^ )], the Friedmann equa- 
tion m) leads to 



'H^ = '^oa'—g* [1 + -^]T 



90 



9* 



(19) 



with 
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is the value of (p at the time of nucleosynthesis. Nucleosynthesis therefore imposes that 
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This constraint, in the framework of quintessence, was first derived by Bartolo and Pietroni pC| . It was also 
shown ||36|1 numerically that in some nonminimally coupled quintessence models the Helium abundance can be 
reduced extending the upper limit on the number of neutrino to 5. Let us also emphasize that very large values 
of |F((/)o)/F(0„„c)| were shown |El| to be consistent with the observed abundances of light elements if (PA/dcjP 
is large enough. Hence, the naive limit (21) can be much more stringent than a detailed (numerical) study may 
show. 



1. Nonminimally coupled scalar field 
As a first example, we consider the case of the nonminimally scalar field introduced in Ref. |l9] for which 



F(0) = l + 2e«/ ((/.), 



(22) 



and we choose the function / to be 



/(</>) = ^^'■ 



(23) 



In Fig. n^, we depict the typical evolution of the equation of state parameter uJq = Pq/ Pq of the quintessence field and 
the time variation of the Newton constant in the minimally coupled and nonminimally coupled cases. As explained 
in the appendices, the equation of state parameter of the quintessence field has two distinct contributions: Wmc that 



appears in the minimally coupled case, and an extra contribution lu^ that appears only in the non minimally coupled 
case, luq, uJmc, and Wj are related by 
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It can be trivially checked that when the tracking solution is reached, 

1 

' 3 



(25) 



Now, a subtlety arises from the fact that depending on the sign of ^, the quantity p/, can be negative. In this case, 
luq does not lie between w^ and Wmc- With our conventions, p^ and ^ are of opposite sign. In both cases one finds, 
however, that the field reaches Planck values as it starts to dominate, i.e., today. Given the form of the coupling 
function (|23|), this means that all the departure from general relativity are mostly felt today. This also implies that 
these departures are stronger in the case of an inverse power-law than in the case of a SUGRA potential, mainly 
because in the latter case the potential is less steep (due to the exponential term) so that the scalar field is more 
stabilized and rolls down slower, implying a smaller time variation of the coupling function F. Another interesting 
feature lies in the fact that coq — ^qL^q changes sign around z ~ 1. 

In Fig. 0, we sum up all the preceding constraints in the plane (fi^,^) both for the inverse power law (m) and 
SUGRA potentials (||). Since ~ Mp, today, F{(J)q) — 1 cx ^ and ^ is then a direct measurement of the deviation from 
general relativity today from which it follows that ^ must be very small today (as first pointed out in Ref . ||2^ ) . The 
constraints are stronger for inverse power law potentials. This is easily understood if one notices that for small ^, 
the equation of state parameter ut is higher in the SUGRA case. This translates into a stronger time dependence of 
the quintessence field, which in turn puts more stringent constraints on the coupling between the scalar field and the 
metric. This remark therefore extends the bounds obtained in Ref. [ p9| . In both cases, the most stringent constraint 
arises from the constraint on the post-Newtonian parameter 7pn and all these constraints become stronger for high 
flq. Another interesting point would be to study how these bounds depend on a. For the Ratra-Peebles potential, 
this was already addressed in Chiba's work: the bound on ^ varies as a~^, which can be checked numerically. For the 
SUGRA potential, the answer is even simpler: the dynamics of the quintessence field varies very weakly with a [ni|. 
Therefore, the bound on ^ are not significantly dependent on a. 
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FIG. 1. Evolution of the quintessence equation of state parameter ujq as a function of the redshift z for Q,q = 0.7 and various 
values of the coupling parameter ^ (left). For negative 5, — | < ojq < ujb in the tracking regime, whereas this condition does not 
hold for positive ^. This is due to the fact that p^ is negative when ^ is positive. For all the models presented here, luq — ojqL ^ 
changes sign around z ~ 1. Evolution of the gravitational constant as a function of the redshift z for various values of the 
coupling parameter ^ (right). The color codes are identical in the two plots. 
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FIG. 2. Summary of the constraints on the parameters (^, fiq) for an inverse power law potential with a = 11 (left) and 
for a SUGRA potential with a = 11 (right). The allowed region lies between the lines of same color. The red lines represent 
the constraint arising from the 7pn post Newtonian parameter [see Eq. IllJ, the green line is from /3pn, the blue line from the 
time- variation of GNefti and the purple line from nucleosynthesis. In both cases, the most stringent constraint comes from 7pn 
and are stronger for high Q,q. 



2. Exponential coupling 
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As a second example, we shall consider the class of scalar-tensor models in which the coupling function A is given 
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with either the inverse power law or SUGRA potentials as defined in Eqs. (|^) and (H). 

To compare with, the standard Damour-Nordtvedt mechanism of attraction of scalar-tensor theories towards general 
relativity, one requires the coupling constant a to drive the field towards its minimum where = 0. In the case of 
quintessence, the field evolves toward infinity at late time and we would need to consider a coupling function a such 
as 



«dn(?) = -Be-''^ 
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that tends to zero at infinity in order to converge toward general relativity. The function A is easily obtained by 
integrating Eqs. (R^, (|27|) to give 
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where the subscript i refers to some initial time. On the contrary, the class of models (26) does not exhibit the double 
attractor mechanism |30| because a and hence, the PN parameter 7pn are constant: one has Oex — —f3, 7pn ^ l~2/3^, 
and /3pJ = 0. therefore, as long as (3 is sufficiently small (we shall take /? = 0.025 in the following), this model can 
be compatible with the Solar system constraints. Note that in this model, the constancy of Oex requires that the 



coupling function F is a polynomial of degree 2, exactly as in the first example. The class of models (27) can indeed 
be easily studied along similar lines. 

In Fig. y|, we show the evolution of the equation of state parameter u; and the Newton constant with time. In 
Fig. H, we sum up all the constraints detailed above in the plane (fi^, B). As already noted, the constraint on the PN 
parameter /3pN is trivially verified, and the constraint on 7pn is also satisfied as long as the parameter f3 in Eq. (|26| ) 
is small. Therefore, only the two other constraints play a role, in contrast with the former case. 
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FIG. 3. Evolution of the quintessence equation of state parameter ojq (left) and of the Newton "constant" Gn^ic (right) as 
a function of the redshift z for ilq = 0.7, f3 — 0.0025, and various values of the coupling parameter B. 
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FIG. 4. Summary of the constraints on the parameters (-B,fi^) in the case of a Ratra-Peebles (left) and SUGRA (right) 
potential with a = 11. The color codes are the same as for Fig. 



B. Supernovae data 



The use of type la supernovae to constraint the cosmological parameters (and hence the claim that our universe 
is accelerating) mostly lies in the fact that we believe that they are standard candles so that we can reconstruct the 
luminosity distance vs redshift relation and compare it with its theoretical value. In a scalar-tensor t heory, we have 
to address both questions, i.e., the determination of the luminosity distance vs redshift relation (Sec. [IB 1 ) and the 
property of standard candle since two supernovae of diff erent r edshift are feeling a different gravitational coupling 
constant and may not be standard candles anymore (Sec. [IB 2). 



1. Luminosity distance in scalar-tensor theories 



To derive the luminosity distance (dL) vs redshift (z) relation needed to interpret the supernovae data, we rewrite 
Eq. (D as 
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where x = 1/(1 + z), and fi", = —K/Hq. The metric of the constant time hypersurfaces is decomposed as 

j,,dx'dx^=dx^+si{x)d'u;, (30) 

where d'^uj is the infinitesimal sohd angle and Sk{x) = [sin(\/i?x)/V^, X? sinh(x\/— -f^)/V— -^] for isT, respectively, 
positive, zero or negative. With these notations, the luminosity distance is given by p^ 
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(31) 



iJ. 5'A'' la in scalar-tensor theories 



The standard lore is to compare this result with type la supernovae data to extract the cosmological parameters 
assuming that they are standard candles in the sense that their light curve does not depend on the supernovae and 
in particular of z. A time varying effective gravitational constant can affect this picture at least in two ways |43| by 
changing 

1. the thermonuclear energy release, since the luminosity at the maximum of the light curve is proportional to the 
mass of nickel synthetized, 

2. the time scale of the supernova explosion and hence the width of its light curve. 



As we pointed out in the Introduction, it was shown in Ref. |£2| that the relation (M) for a flat Cold Dark Matter 
model with a cosmological constant (ACDM) in the framework of general relativity can be mimicked by a scalar-tensor 
theory with y = A = up to a given redshift. But, this study did not include the fact that in scalar-tensor theories 
one, as will be shown here, can not directly use the luminosity distance vs redshift relation inferred in the framework 
of general relativity. Hence, when comparing any scalar-tensor theory to supernovae type la (SN la) data, one needs 
a modified magnitude vs redshift relation taking into account the effects listed above. The goal of this section is not 
to explain the recent SN la data by replacing the cosmological constant by a scalar-tensor theory but rather to know 
how to deal with these data in such a framework and to evaluate the effect of the coupling on the precision of the 
determination of the cosmological parameters. 

To discuss these issues, we recall a simple "one zone" toy model Q,^^ of an expanding sphere of uniform tempera- 
ture T and density and of radius and mass i?e,„ and M^^^ for the supernovae light curve which encapsulates the main 
features of dependence in the gravitational constant G, even if this model is nothing but a toy model. The observed 
light curve is obtained from the nonadiabatic evolution of the thermal energy £"11, — AimT^R^^^/S of the radiation 
dominated envelope 

£;„, + 47rPi?i^i?,„, = -L + L^, (32) 

where L ~ -Eth/'^ditt is the bolometric luminosity. L^ = L*oCxp[— t/r*] is the radioactive energy input (from 5gNi -^ 
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Fe, where r* represents the corresponding nuclear process half life time). P = aT /3 is the radiation 



pressure and the thermal energy E'th and the photon diffusion time scale Tam is given by 



Tditt — 3- 



.Rt 



X. 



47ri?,„, 



(33) 



where K^h is the Thomson opacity and A.^ the photons free mean path. Since the temperature scales as l/i?env, it 
follows that 



^th — 



Mq Rq 

2 i?c.„ 



(34) 



where Rq and AIq are respectively the radius of and the mass of the progenitor. Equation (|3^) yields the equation for 
the luminosity L 



^L=i-L + L..)^, (35) 

^ "-onv ' diff 

where t^.jj. is the initial diffusion time scale when /J^nv — Ro- If the envelope expands at a constant velocity u^^p = i?e„v, 
an analytic solution to Eq. (|5|) was found by Arnett |Q as 

L = Ue-''^'='> +L°n{x,y,w), (36) 



where, from Eqs. (|33|) and (34), Lp ~ -Eth/'^'ditt — 27r_Ro/9'*Th is the luminosity at the maximum of the light curve, 



u{x) = wx + x'^, X = t/r, y = t/2t^, u) = t/t^^^^, and r is the characteristic time of the SN given by 

r^ = 2f^rL, (37) 

and the function f2 takes the form 

J7(a;,y,w) = e-"("W {w + 2z)e~^'y+'''^'^ Az . (38) 
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The expansion velocity can be obtained via the conservation of energy as 



v^ =2^. (39) 

Mo ^ ' 

In order the model to be consistent, we are therefore obliged to assume that M^-, < ^Mq. From the above equation, 
it follows that r behaves as 

r^.^fMV^ (40) 

Assuming that the progenitor has the Chandrasekar mass Afch, -Ro is the Chandrasekar radius, and assuming that 
the mass of nickel scales as Mq , we deduce that 

L....-LZt[^) , r = r-(^^j , (41) 

since the total energy release is proportional to mass of nickel formed which is assumed to scale as the progenitor 
mass and thus as G""^/^. 

Such a toy model can, at that stage, describe both SN I and SN II. But, for SNIa, the progenitor is a white dwarf 
and it follows that i?o ^ 5 x 10"^ km. Then, since 

and 

r~71— = ( , "'''' , ] , (43) 

Mo VlO^kms-i/ ' ^ ^ 

we can conclude that Lp ~ and that lu ^ 10~^. We then choose the typical M^ value of the parameters for SN la 
to be 

Lp = 0, w = 0. (44) 

It follows that the luminosity curve is well approximated by L = L'^n{x, y, 0) which can be integrated analytically to 
give 

L = L2 (e'^'"'" ~ e^' + i\/^ye-(^'+^') [Erf(ia; - iy) + ETi{iy)]\ , (45) 

where Erf is the error function and where the nuclear rate are given by 

10 



cobalt: L, = 1.5 x 10'*^MNiergs"\ r, = 111 days, 
nickel: L, = 8 x lO'^^MNjergs^^ r, = 10.1 days. 



(46) 

(47) 



In Fig. g, we depict a standard light curve obtained with this one zone model and the light curves when G is 
increased respectively by 10% and 20%. This light curves are compared to the ones obtained by a variation of the 
nickel mass synthetised. The decaying branch is mainly sensitive to the mass of nickel so that an increase of the 
gravitational constant implies that the light curve has a lower maximum and then tends asymptotically towards the 
light curve of a supernovae with lower nickel mass. As a conclusion, if G^it is 10% larger than G all other parameters 
being unchanged, then the luminosity at maximum will be slower by about 15% and the time scale r will be smaller 
by 7.5% and the light curve will be narrower. As was realized, supernovae are not exactly standard candles but, 
thanks to the correlation between the time scale of the light curve and the peak luminosity (larger curves are brighter 
while narrower are fainter), the dispersion of 0.15-0.2 mag can be corrected by the use of a stretch factor H]. 

A variation of G affects both the amplitude of the peak and the time scale and, for instance, makes it narrower and 
fainter if G grows in the past. Again, they can be calibrated to extract the luminosity distance since, once the model 
is specified, the dependence of the correction due to the variation of G is known. The magnitude vs redshift relation 
then takes the form 



m(z) = TWo + 5 logd^iz; f^o, f7° , C] 



15 , G,„ 



(48) 



if we just take into account the effect on the peak luminosity (assuming that a stretch factor has been applied 
yet). Our result is compatible with the relation obtained in the case of a Brans-Dicke theory p3| for which G = 
(4 + 2a;)/(3 + 2ijj)(j)~^. Such a correction was also argued in Ref. [Q where it was assumed that the peak luminosity 
scales as GZa and that G^^ — G(l + tHo)"^ leading to what was referred to as a "G correction" in the magnitude 
vs redshift relation of the form Amc — 2.5m7log[l + t{z)], t{z) being the look-back time. Under these hypothesis 
these authors showed that an increase in G in the past can reconcile the SN la data with an open, A = cosmology. 
But, what was not shown is t hat such a variation can be cast into a scalar-tensor framework while respecting all the 
constraints described in Sec. II A and particularly respecting the positivity of the energy of the graviton [22 



The 

phenomenological exponent 7 is not a free parameter and has to be determined from the theory. Here, we claim, on 
the basis of our toy model, that 7 ~ 3/2. 

Equation (|4^) gives the general magnitude vs redshift relation in scalar-tensor theory and has, to be self-consistent, 
to be used when comparing an extended quintessence model to the data. Indeed, the one zone toy model presented 
above may be thought to be very rough but is qualitatively correct Q . We have to stress that we neglected the effect 
of the scalar field in the SN la dynamics and assumed that its only effect (or at least dominant) was to change the 
value of the gravitational constant. We also point out that we used Getf instead of GN„ff in Eq.(H^) but, as shown in 
Ref. WTl they may not differ from more than 10%, which is a 1% effect in our magnitude vs redshift relation. 
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FIG. 5. Variation of the light curve of SN la with respect to an increase of the gravitational constant in the simple "one 
zone" toy model presented in this article. A stronger gravitational constant accounts for a fainter and narrower light curve (see 
Eq. pil)). The rising part is mainly sensitive to the progenitor mass whereas the decaying part is mainly dictated by the mass 
of nickel. Hence each light curve computed with a higher Newton constant tends asymptotically toward the light curve of a 
supernovae with smaller nickel mass but with standard value of the Newton constant (left), and vice versa (right). [Note that 
in this last case, we have to ensure that the total nickel mass does not exceed half of the total mass in order for Eq. (p9|) to be 
valid.] 

3. Application to two examples 

In Fig. |6| we show the luminosity distance vs redshift relation for the two models presented above under the 
assumption that the universe is spatially flat (i.e., rj^^t+r^Q = 1). In the case of the nonminimally coupled quintessence 
field, it can be noted that the deviation is more important for positive ^ than for negative ^. Note that in any case, 
the net effect of the coupling of the luminosity distance is small (the most extreme values of S, and B plotted here 
are already ruled out). However, when fitting the supernovae data, one must also take into account the change in the 
SN absolute luminosity, which can be as large as 30% for the largest variation of G allowed by nucleosynthesis. We 
also note that the distance vs redshift relation is more sensitive to the coupling S, for smaller values of the exponent 
a. From the Fricdmann equation (0) and Eq. (0), we have 



Arj oc F. 



(49) 



Moreover, for constant 0, increasing ^ increases -F, and therefore the distance traveled by a photon is increased. This 
explains why the quantity cJl — d^ increases with ^ in Fig. g. Unfortunately, the same reasoning does not apply for 
the exponential coupling because the hypothesis that the value of the quintessence field (p does not significantly vary 
at fixed redshift when one varies the parameter B is no longer valid. (Should the field be at rest today, then the 
hypothesis would still be vahd, but this is not the case, as shown in Figs. || and||.) 

In Fig. ^ we present the distance modulus vs redshift relation taking into account the supernovae luminosity 
correction due to a variation of the Newton constant. We compare the amplitude of this correction to the standard 
ACDM model. It appears that the allowed values for ^ do not lead to a significant modification of the distance 
modulus. This is related to the fact that this model is very strongly constrained by the solar system data. On the 
contrary, the exponential coupling model allows for a larger deviation of the distance modulus, mainly because it 
is not limited by the post Newtonian constraints. More generally, any model tuned so to evade the solar system 
constraints can lead to large deviation in the magnitude vs redshift relation. 
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FIG. 6. Luminosity distance vs redshift relation for different value of ^ (left) and B (right). Both are compared to the 
luminosity distance vs redshift relation d* for a matter dominated, flat universe. As noted previously, the deviation is more 
important for positive ^ than for negative ^. 
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FIG. 7. Distance modulus vs redshift relation for different models. A comparison between a Standard Cold Dark Matter 
(SCDM) model, a ACDM model, and a Quintessence Cold Dark Matter (QCDM) model are shown. Also shown are the 
maximum allowed deviation from a standard QCDM model when considering the two models described in the text. The 
exponential model allows for large deviations (roughly equivalent to those obtained by adding SQq ~ 0.05 in a minimally 
coupled QCDM model), whereas the non minimal coupling is much more constrained. The right panel represents a zoom of 
the left panel for high redshift. 



C. Angular distance 

Another important observational quantity is the angular distance, cIa, which, with the notation of the previous 
paragraph, is given by ' 



d,= 



[i + zy 



(50) 



This distance relates the size of an object to the angle under which it is observed. Naively, the variation of d^ can be 
related to the position of the first cosmic microwave background acoustic peak. Note also that (^^{z) also enters the 
probability of lensing y] . 

In Fig. pi we depict the variation of dj^{z) both in the nonminimally coupled scalar field models and in the exponential 
coupling case. Indeed, a complete study of the CMB anisotropies will be presented in the following section and these 
results are just a hint of how the position of the first acoustic peak will be affected in these models. The relative 
shift in the position of the acoustic peaks is simply given by the asymptotic part of the curves of Fig. g. As for the 
luminosity distance, the effect of the coupling is rather small. However, the variation of G will lead to an important 
modification of the peak position as we shall see later. 
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FIG. 8. Angular distance vs redshift relation for different values of ^ (left) or B (right) compared to dX, the angular distance 
vs redshift relation for a matter dominated, flat universe. The relative shift in the position of the acoustic peaks is simply given 
by the asymptotic part of the curves. We assumed that SIq — 0.7 and used the SUGRA potential with a — 11. 

III. PROPERTIES OF THE PERTURBATIONS 

The background and perturbation equations in a scalar-tensor quintessence model are presented in the appendices. 
Not that these equations are not restricted to the case of the nonniinimally coupled scalar field and can be extended 
easily to the exponential coupling example by simply setting 

since, as emphasized in the Introduction the splitting (||) is completely general as long as the function / has not be 
chosen. Technically, one obtains a set of equations slightly more complicated than in the minimally coupled case, but 
which can still easily be solved numerically. We shall first discuss the main properties of the scalar field perturbations 
before turning to cosmic microwave background properties in the next section. 

As for the minimally coupled case, one of the main concern is about the choice of initial conditions for the pertur- 
bations. It happens that this problem is not relevant since the long wavelength scalar field perturbations follow an 
attractor for the following reason. First, let us consider the unperturbed and perturbed Klein-Gordon equations 



AScp + 2rL5<i) + a^V"5<i) - 2X0 - 



2n^ + aV = -^a^Rf, (52) 

-£,{a^Rf"5(t) + a^6Rf'), (53) 



A— -X-4HX 
tL 



where $ is one the two Bardeen potential and X is a combination of the two Bardecn potentials (see the appendices 
for more details). For an inverse power law potential, one can show that there exists some particular solution to 
the unperturbed Klein-Gordon equation, the so-called tracking solution |j], for which the field evolves according 
to some power law (the exponent of which depends on the exponent of the potential). Then, the stability of this 
particular solution — which is the most useful feature of the quintessence scenarios — is determined by the properties 
of the "perturbed" expression of the Klein-Gordon equation. This "perturbed" equation actually describes the small 
departures of an homogeneous field from its tracking solution, and reads 

5(j) + 2m(j) -t- a^V"5(l) = -^a^Rf'Scj). (54) 

Obviously, this equation is extremely similar to the above "real" perturbed Klein-Gordon equation: we simply have 
neglected the metric perturbation terms, and considered the large wavelength limit fc ^- 0. The very nice feature of 
inverse power law potential is that the solutions to this equation tend to 0. Therefore, the solutions of the homogeneous 
part of the real perturbed Klein-Gordon equation in the large wavelength limit tend to 0. Now, when one takes into 
account the metric perturbations, remembering the fact that they are constant in the long wavelength limit (see, e.g. 



Ref. 1 49 1), it appears that the relevant quantities, such as Sm_c (the density contrast of the quintessence fluid) will 
tend toward constants which will simply be linear combinations of the metric perturbations. Although not explicitly 
stated, this is what was shown in Ref. [ p7[ . Finally, in the short wavelength limit, the Laplacian term ensure that the 
field follows a wave equation, which is actually damped by the expansion. 

As a conclusion, as long as the field is subdominant (which is the case when one fixes the initial conditions and 
in the tracking regimes), its perturbations follow an attractor, thus solving the problem of the choice of the initial 
conditions, and the field does not show any unstable behavior, as illustrated in Fig. 0. 
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FIG. 9. Evolution of the scalar field density contrast as a function of redshift for a non minimally quintessence model 
(SUGRA potential with a = 11, power law coupling, and ^ = 0.06). We have represented the evolution of the Fourier mode 
of the quantity Sf^Q{k,z) for the wave number k = 0.1/iMpc~^. By comparison, we have also shown the Bardeen potential 
$. The quintessence field reaches the attractor around z ~ 10^^. After this epoch, the quintessence density contrast reaches 
a value depending only of the potentials $ and <&. It subsequently evolves only when the mode has entered into the Hubble 
radius, where it experiences damped oscillations (after z ~ 10^) around another value also depending on $ and ^. The 
quintessence density contrast can be very large before having reached the attractor. This is, however, not a problem, because 
the quintessence density contrast JImc is small at this epoch. Note that the Bardeen potential has a nontrivial behavior at 
late times both because of the domination of the quintessence field (which explains the decay after z; ~ 3) and because of the 
variation of the Newton constant. 

IV. COSMIC MICROWAVE BACKGROUND AND MATTER POWER SPECTRUM 

Appendix gives the derivation of the modified equations for the evolution of the cosmological perturbations. 
As explained in this appendix, the two most important modifications induced by the quintessence field are (i) a 
modification of the background equations which lead to the domination of the scalar field today and (ii) a subse- 
quent modification of the late time evolution of the gravitational potentials. We have implemented the scalar-tensor 
quintessence equations in a Boltzmann code (which uses the line-of-sight integration method, see, e.g., Refs. [pO| , pl[ ) 
and computed matter power spectra as well as CMB anisotropics. We shall now discuss the observable consequences 
of these models. 

Two effects of a quintessence (minimally coupled or not) have already widely been discussed in the literature |^,^ . 
These are the following: 

• a modification of the angular scale of the peak structure: By modifying the expansion rate of the Universe at 
low redshift, the quintessence modifies the usual angular distance vs redshift relation (see Fig. |^). This induces 
a global shift of the acoustic peak structure. Of course, if one adds a quintessence field while keeping fixed the 
Hubble parameter today, then this is equivalent to modifying the redshift of equivalence between matter and 
radiation, which of course has some influence of the acoustic peak structure (the same problem occurs for a 
cosmological constant). From Fig. ||, we would have expected the first acoustic peak to be shifted to smaller 
multipoles since the diameter angular distance is smaller when ^ > 0. This is indeed not what is observed on 
Fig. O. This is because another, more important effect is that a variation of G modifies the Friedmann equation 
at early times, and therefore the age of the Universe (and, hence, the sound horizon) is modified. For example, 
when 1^ > 0, G is larger at earlier time (see Fig. nh, and the age of the Universe is smaller at recombination, 
which shifts the peak structure towards higher angular scales. 

• a boost of the Sachs- Wolfe plateau: When the field starts to dominate, the density parameter of the ordinary 
matter decay, as well as the gravitational potential. This leads to the possibility of energy exchanges between 
photons and gravitational field, hence producing the so-called Integrated Sachs- Wolfe (ISW) effect which will 
boost the anisotropy spectrum of scales larger than the Hubble radius at the epoch of transition between matter 
and quintessence. The cross correlation between the Sachs- Wolfe and the ISW terms are difficult to compute, 
and, as long as the ISW term is not too important, one can either have a higher or lower first peak. 
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For any realistic model, the field is by far subdominant at recombination. As a consequence, the acoustic peak 
structure of the CMB anisotropies will not directly be affected by the field dynamic. However, two new effects arise 
when one introduces a non minimally coupled quintessence field: 

• a modification of the amplitude of the Silk damping: At small scales, viscosity and heat conduction in the 
photon-baryon fluid produce a damping of the photon perturbations [ pl| . The damping scale is determined 
by the photon diffusion length at recombination, and therefore depends on the size of the horizon at this epoch, 
and hence, depends on any variation of the Newton constant throughout the history of the Universe. 

• a modification of the thickness of the last scattering surface: In the same vein, the duration of recombination is 
modified by a variation of the Newton constant as the expansion rate is different. It is well known that CMB 
anisotropies are affected on small scales because the last scattering "surface" has a finite thickness. The net 
effect is to introduce an extra, roughly exponential, damping term, with the cutoff length being determined by 
the thickness of the last scattering surface ||5^ . This thickness is of course determined by the duration of the 
recombination process. This process is essentially fixed in redshift units, as it is dominated by Boltzmann factors 
(and hence, by photon temperature) and not by standard cosmological parameters. However, when translating 
redshift into duration (or length) , one has to use the Friedmann equations, which are affected by a variation of 
the Newton constant. The relevant quantity to consider is the so-called visibility function g, defined as 



giv) 



(55) 



with f being the differential opacity, and rjtoday = 0. In the limit of an infinitely thin last scattering surface, 
r goes from oo to at recombination epoch. For standard cosmology, it drops from a large value to a much 
smaller one, and hence, the visibility function still exhibits a peak, but is much broader. A few examples are 
given in Fig. 10 for several values of the coupling parameter ^. As it can be seen, the height of the visibility 



function varies with ^, and therefore, its thickness also varies as the integral of this function is 1. 




FIG. 10. Modification to the visibility function due to non minimal coupling between a quintessence field and curvature. 
Minimally coupled models are identical to ACDM models since the expansion history of both model is strictly identical at high 
redshift, whereas important differences can occurs when one introduces a non minimal coupling. For a better readability the 
various visibility functions have been slightly shifted upward. 

• a modification of the sound horizon at recombination : as for the photon diffusion length, the sound horizon 
depends on the size of the horizon, at can significantly change if the Newton constant is different from now at 
this epoch. The net effect is to produce a shift in the Doppler peak structure. 

More quantitative estimates of all these four effects are difficult to compute analytically, since there is absolutely 
no reason that there exists any simple solutions to relevant equations. However, as already noted for the supernovae 
luminosity distance, the effect of the non minimal coupling one the angular size of the acoustic peaks is small, except 
for unacceptable large value of the parameters. On the contrary, the ISW effect exhibits a much stronger dependence 
on ^ and B, as shown on Fig. |ll|. One also notices significant modifications in the damping at high multipoles. 

On Fig. n3, we have plotted the corresponding matter power spectra. Two important observable effect arise here. 
First, the normalization of the spectra changes because the usual Sachs- Wolfe formula ST jT = <i>/3 which relates the 
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CMB anisotropy to the matter power spectrum amplitudes does no longer hold because of the ISW effect. This is 
particularly obvious for large, positive values of ^. Second, the maximum of the power spectrum is shifted by the 
coupling. This is because this maximum give the scale that enters into the Hubble radius at the epoch of equality 
between matter and radiation. In all the cases presented here, this epoch always corresponds to the same redshift but 
not to the same cosmic time because the Friedmann equations are modified. Therefore, the relation H{z) is different 
in all these models. 

In practice, the power law coupling ( |23| ) is already very strongly constrained by the PN parameters. For all 
acceptable values of the parameters, the CMB anisotropics show a negligible deviation from the minimally coupled 
case (note that on Fig. 12, the plotted non zero values of ^ are one order of magnitude larger than allowed values inferred 
from Fig. g). On the contrary, for the exponential coupling where the PN constraints are fulfilled by construction, 
the CMB anisotropics can play a significant role in constraining (or measuring) the coupling parameters. 
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FIG. 11. Cosmic microwave background anisotropies for the SUGRA potential (q = 11) with a nonminimally coupled scalar 
field (left) and with an exponential coupling (right), respectively, for different values of ^ and B. We have assume a locally flat 
Universe. The strongest features are the boost of the decay of the amplitude on large scales due to the Integrated Sachs- Wolfe 
term, as well as a slight shift of the acoustic peaks structure to higher multipoles. 
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FIG. 12. Matter power spectrum for the SUGRA potential (a = 11) with a nonminimally coupled scalar field (left) and 

with an exponential coupling (right), respectively, for different values of S, and B. We have assumed a locally flat Universe. 

The unusual features are the variation in the normalization of the power spectra once COBE-normalized, and the shift in the 

maximum of the spectrum whereas the epoch of equality is the same in all these models. 
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V. CONCLUSION 

In this article, we have investigated the interpretation of cosmological and astrophysical observations in the context 
of scalar-tensor quintessence. In this class of models, the quintessence field induces a time variation of the gravitational 
constant and the post-Newtonian constraints and the constraints on the time variation of the constants of nature where 
first imposed, as well as the restriction arising from nucleosynthesis. 

We then focused on the supernovae and cosmic microwave background dataset. Concerning the supernovae, we 
extracted, using a simple toy model, the dependence of the maximum of the light curve on the gravitational constant. 
We also explained the various effects of this coupling (mainly the modifications it induces of the Newton constant) on 
CMB anisotropics and structure formation. All these features on the CMB angular power spectrum and the matter 
power spectrum are physically well understood. 

All these observational constraints were applied all along the article to two class of models: (i) a nonminimally 
coupled quintessence field and (ii) and exponential coupling. We showed that the first class of models was very 
constrained mainly because the deviation from the general relativity was fixed and that the constraints were more 
severe in the case of an inverse power law potential than for the SUGRA potential. The second class seems to leave 
more freedom for parameters and could possibly be constrained using the CMB. 
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APPENDIX A: LAGRANGIANS 

The general action (ph with the ansatz (ph can be decomposed as 

S= [[€, + £mc + C^. + C^J /Md^^, (Al) 

representing respectively the Lagrangian contribution for the gravity (g), the minimally coupled scalar (MC) field, 
the coupling (co), and the matter (mat) and given explicitly by 

A = ~^R, (A2) 

Cmc = ^g'^'D^cf^D,^ - 1/(0), (A3) 

/:co = -ei?/(0), (A4) 

where Da is the covariant derivative of the metric gafs- By varying the Lagrangian £q, we obtain the following 
stress-energy tensor 

T^^ = Da^Dpcf> ~ ^D^<l>D^(l>gap + V{cf>)gaf3. (A5) 

Varying the Lagrangian C^^ describing the coupling, one obtains a stress-energy tensor that can be separated into two 
parts 

r~^ = ~2^Gapf, (A6) 

n^^-^^gapD-DaDp)/, (A7) 

where Gaf3 is the Einstein tensor and D = Dfj^D^^. We have separated the coupling term into two components to single 
out the part (labeled k) which can be absorbed in a redefinition of the Newton constant. It follows that the Einstein 
equation takes the general form 
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G^p = n Y. ^i/3 (A8) 



/ = MC,£,K,mat 



or, equivalently, with the use of Eq. (g) 



Go.0 = n^u Y. T^p, (A9) 



/ = MC,5,mat 



with 



K..ti = — = . (AlO) 



APPENDIX B: BACKGROUND SPACETIME EQUATIONS 

By varying the Lagrangians with respect to the scalar field, we obtain the Klein-Gordon equation 

□0 + F' = ~^Rf (Bl) 

or 

'(j) + 2W0 + a^V = -^a^Rf, (B2) 

where R stands for the scalar curvature and a prime denotes a derivative with respect to (j). Alternatively, the usual 
conservation equation (-D^jT"^" = 0) for the matter has to be completed by the conservation equation of each of the 
stress-energy tensors defined above 



D^TZ = -2C[-]<?^"i^^/ ^ QSc: (B3) 



,R, 

D^T^^'^ = -2^[-Rt^'-]D^f = Q^, (B4) 

D^rr = -2^[GnD^J = Q'i (B5) 



where the Q" can be interpreted as a force term |Q. It is straightforward to check that, as expected, the sum of the 
right-hand sides of these thre e equations vanishes. 

When developing Eq. (A9) with the metric (^, we obtain the Friedmann equations 



2,{re + K) = K^,,a' Y Pf' (B6) 

/ = mat,MC,,^ 

2{n^ + K-n) = -K^,,a^ Y (Pf+Pf), (B7) 

where we have introduced the quantities 



/ = mat,MC,,^ 



a^Mc = ^0' + aV, (B8) 

a^P,,c = ^0' - aV, (B9) 

a^p, = -em, (BIO) 

a^P, = 2(if + nf). (Bll) 

Pmc and Pmc correspond to the minimally coupled part of the scalar field energy density and pressure, p^ and P^ to 
the nonminimally coupled part which can not be absorbed in K^ff. Note that these two equations are not completely 
straightforward to solve since there are terms proportional to Ti, K, Ti. in p^, P^. For example, defining 

(B12) 
(B13) 
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n,^ 


-K + 


Kcfffl^ 


>: 




3 


Pf 




V 


/ 


= mat,MC 




£5 = 











one gets 



Ti = TCf 



VT" 



from which one can compute p^. For Pj, one can use the foUowing formulas: 

/3 



p = p. 



1 + 3/3 



Pe 



3P,+ Y. (Pf'^Pf) 



/ = mat,MC 



We also set 



Pq =Pmc +P5, 



(B14) 

(B15) 
(B16) 

(B17) 



(B18) 
(B19) 



which are the contribution of the scalar field entering the right-hand side of the Friedmann equations, and 

a^p^ = -6$ {H^ + K) f, (B20) 

a^P, = 2^ hn + n^ + k\ /, (B21) 



which are the quantities entering the redefinition of k. With these notations, the conservation equations (B3) 
take the form 



Pmc + 37Y (p^ 






n^ + K + n 






P. + sn [p^ 



Pj = ^ [n' 



K 



(B22) 
(B23) 
(B24) 



and it can be checked, as expected, that the sum of the right-hand sides of Eqs. (B22)-(B24) vanishes. 

In practice, it is not difficult to solve this set of equations [apart from the subtleties described in Eqs. ( |B12| )- 
(B17)]. The main problem arises from the fact that we know a priori neither the values of the energy scale of the 
potential (l|),(0), nor the bare Einstein constant k. All what we know (or impose) are the value of fig and Ketf today. 
We are therefore obliged to use standard routines Q to find a solution that converges towards the desired values 
for riq and k^h today. In order to do so, one is obliged to provide a reasonably good starting point for k and M. In 
practice, one can consider 






(B25) 
(B26) 



which usually converges for reasonable values of the parameters. When this starting point fails to converge, we are 
usually in a region of the parameter space already ruled out by the astrophysical constraints detailed in Sec. [I A. 



APPENDIX C: PERTURBATION EQUATIONS 

We consider the general form of a perturbed cosmological spacetime 



ds' 



(1 -I- 2A)dr]^ + 2{ViB + B,)d7]dx' - (7^ + 2C7y -t- 2V^VjE + 2V(,£:,) -f 2Eij)dx'dx^ 



(CI) 



where 7^ is the metric of the homogeneous spatial sections and V^ its covariant derivative. The vector quantities Bi 

and Ei are divergenceless {S/^Bi = VEi = 0), and the tensor Eij is divergenceless and traceless. 
For any fluid, the stress-energy tensor can be decomposed as 
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Tap = (-P + p)UaU/3 - PgaP + ^ap, (C2) 

where liap is a traceless tensor orthogonal the the four- velocity m". For a perfect fluid, we have IIq^ = 0. We then 
define 

-a^{P + p)v, = {P + p)uogijSu^ + SUo,, (C3) 

a^Pn,,=6Il,,, (C4) 

where Vi and TTy have to be seen as three-dimensional quantities whose indices are raised and lowered by the metric 
jij. Using these quantities, the stress-energy tensor perturbation takes the form 

SToo = pa^ (S + 2A) , (C5) 

JTo, = pa^ [V,B + B, - (1 + Lu){V,v + v,)] , (C6) 



dT,j = Pa 



SP _ = 

hij + —f^lij + AyTT + D(^i1Tj) + TTij 



(C7) 



where hij = 2C^ij + 2ViV jE + 2V(i-Ei) + 2i?y, A^ = ViVj — 7^ A/3, A = ViV* being the Laplacian, and tt, tt^ and 
Tfy are the scalar, vector and tensor components of the anisotropic stress tensor. 

We work in the gauge in which B = E = and Ei = and introduce the gauge invariant perturbation variables, 
labeled with a superscript (j, 

5i=^-P^P^B + E), (C8) 

P P 

v^ = v + E, (C9) 

vl = v,+t,, (CfO) 

and the four gauge invariant gravitational potentials 

^ = A + HiB + E)-{B + E), (Cll) 

<P = -C + niB + E), (C12) 

V,=B,+%. (C14) 

(Note that X can be expressed in terms of $ and ^ and their time derivatives. It will, however, prove useful to work 
with these three quantities.) The pressure perturbations are related to the density contrasts by 

^ = ^<5«+r, (C15) 

where cj. = P/ p is the (adiabatic) sound speed and F is the entropy perturbation. As we shall see later, F^ and F^c 
do not vanish. We also introduce the convenient flat-slicing gauge in which C = E = and Ei = (and thus where 
A = X and B — ^/7i) and where the gauge invariant density contrast and pressure perturbations, labeled with a 
superscript b, are given by 

^^.^-^^=^« + ^* (C16) 

P P n p n 

SP^^SP^ P^ 

P P PH ^ ' 



Note that in Eq. (pl6) we do not use the conservation equation to express p in order for this definition to be valid 



also for coupled fluids. The velocity perturbations are identical in both gauges. 

The scalar field is decomposed as (j){r]) + S(j) and we introduce the two gauge invariant variables, respectively, in 
Newtonian and flat-slicing gauge 

S(f>* = S(j)~'<j){B + E), (C18) 

Sd,' = Scf,^ + ^^. (C19) 
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It is indeed not the purpose of this Appendix to rederive all the equations of perturbation; our aim is to take into 
account the nonminimally coupled scalar field. For that purpose, we first compute Klein-Gordon equation. Then, we 
compute the perturbed stress-energy tensor to deduce the Einstein equations, using the standard scalar-vector-tensor 
decomposition. 



1. Perturbed Klein-Gordon equation 

The perturbed Klein-Gordon equation reduces to 



„■ b 



$ 



A— - X - 4nx 
n 



AScf)" + 2Hd(p' + a^V"5(f}' - 2X0 + </- 

where we have used the scalar curvature R 

a^R = -Q{n'^ + K + n), 
and of its perturbation dw , in the flat-slicing gauge 

a^5R^ = 2A(vE' - 2$) + 6(27i:^X + 27i:X + HX) 



= -i{a'Rf"54f +a^5R'f'), 



(C20) 

(C21) 
(C22) 



2. Perturbed fluid quantities 

After some manipulations, it appears that it is possible to arrange the stress-energy tensor T^S so that a minimally 
coupled scalar field essentially behaves similar to a fluid: 



Pmc = ]^D^,<j)D^'c^ + V, 



D"6 



^d;jd^ 






From these expressions, it is then possible to compute density and pressure perturbations, etc. 

a^5pl,^ = 0(5/ + a^V'Scf}' + ^^X, 



a'-^P^c 



a^{Pmc + Pmc)vl 






a^iPM 



= i^V, 



lie 1 PmcJ^mc i 

a PmcT^ij = U. 

Similar expressions can be found for the ^ part. For the density and the pressure, we have 



Pi 



i(-6en/ + 3V7), 



with 



n- J(6en./ + V7), 



y=^ [47^.^^^ - {g,.Tn^] 



^if-2'Hf) 



(C23) 
(C24) 
(C25) 
(C26) 

(C27) 

(C28) 

(C29) 

(C30) 
(C31) 

(C32) 
(C33) 

(C34) 



As for the minimally coupled part, we set 
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from which we deduce an anisotropic stress H" : 



jja/3 ^ jj.. 



T-Q/3 



r"-7ff 



4- ^« 



5m.] - VF (^^i^f - i^"' 



One can easily check that H"'' = and 5Wf = 0. Using Eqs. (C3),(p|),(C6),(C7), we finally obtain 



a^5p\ = 2^ 



a^SPl = 2£_ 



a\P,+p,)v\^2£, 



x^ 



-snsf +6nfx + AUf-f^ 

6f + mf - ^A Uf - j^ -f±- 2/X - 2H/X 



A. $ 



Sf' + n5f + if-2nf)- + fx 



n 



a^P.n^ = -2i (f% - Sf 



H 



a'P.nl = -2C/y„ 
a'P,% = ~2tfE,,. 



Finally, the perturbation of T^o gives 
In the next section, we shall inject these results into the Einstein equations. 



(C35) 

(C36) 

(C37) 

(038) 

(039) 
(O40) 
(041) 

(042) 
(043) 

(044) 



3. Einstein Equations 

With these quantities and the standard expressions of the perturbed Einstein tensor, which can be found, e.g., in 
Ref. [p5|, we can give the perturbed Einstein equations 



/ = mat,/t,S,MC 

in the scalar-vector-tensor decomposition. 



/ = mat,Kji,!;,MC 



(045) 



a. Scalar modes 



After some long but straightforward manipulations, we finally obtain 



2(A$ - m^X) = m^ ^ ^fis) 



/ = mat,£,MC 



Sk!' 



'H^X+ iH-H^]<^ 



m^ Y, nf{i + Lo,)nv), 



/ = mat,5,MC 



* - * = 37^2 J2 Qfiv.TTf, 



/ = mat,,c^MC 



n^X + 2'HX + HX+ - A(\E' - $) 



2 V- r. ^P\ 5k\., 



/ = mat,,c^MC 



(046) 

(047) 
(048) 

(049) 
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We therefore end with a hnear system relating the four metric perturbations $,^$,X, X to the fluid perturbed 
quantities, some of which depending on the metric perturbations, [cf Eqs. ( |C27| )-( p43| )]. This system can easily be 
solved numerically in order to extract the metric perturbations, which in turn source the perturbed mass conservation, 
and Euler and Klein-Gordon equations. Apart from the different expressions for the metric perturbations, the other 
fluid conservation equation are not modified by the presence of the quintessence field. Hence it is quite straightforward 
to implement this in any preexisting numerical code. 



b. Vector modes 



In the same way, we obtain 



1 



^{A + 2K) + 2{n^ + K - n) 



2nWi + nf^ ^ 'SH^ Yl ^i^M' 



(C50) 
(C51) 



f=mat,^,MC 



which is then easily solved. 



c. Tensor modes 



The equation of evolution of the gravitational waves reads 

E,j + 2HE^j + {2K - A)E^j = SH^ ^ ^/W/^ 



ij- 



(C52) 



/ = mat,^,MC 



The presence of the term proportional to Eij in Eq. ( p43| ) indicates that the damping rate of the gravitational waves 
is different. This result is identical to the one we had obtained in an earlier work [ |l3| . 
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